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Abstract 

We have used 87 ptT 1 of data collected with the Collider Detector at Fermilab to 
search for new particles decaying to bb. We present model-independent upper limits 
on the cross section for narrow resonances which exclude the color-octet technirho 
in the mass interval 350 < M < 440 GeV/c 2 . In addition, we exclude topgluons, 
predicted in models of top color-assisted technicolor, of width T = 0.3M in the mass 
range 280 < M < 670 GeV/c 2 , of width T = 0.5M in the mass range 340 < M < 640 
GeV/c 2 , and of width T = 0.7M in the mass range 375 < M < 560 GeV/c 2 . 
PACS numbers: 13.85.Rm, 12.38.Qk, 14.70.Pw, 14.80.-j 
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In this paper we report on a search for new particles decaying to bb. In addition 
to a model independent search for narrow resonances, we perform a specific search 
for resonances from top color-assisted technicolor |]J . Since electroweak symmetry 
breaking is associated with the origin of fermion masses, the large mass of the top 
quark suggests that the third generation could contain clues about the origin of elec- 
troweak symmetry breaking. Top color-assisted technicolor is a model in which the 
top quark is heavy because of a new interaction. Topcolor replaces the SU(3)c of 
QCD with SU(3)i which couples only to the first two quark generations and SU(3) 2 
which couples only to the third generation. The symmetry S77(3)i x S77(3)2 is bro- 
ken, resulting in the familiar SU(3)c of QCD and an additional SU(3) which couples 
mainly to the third generation. In addition, there is a U(l) symmetry added to keep 
the b quark light. The additional SU(3) symmetry gives rise to a color-octet gauge 
boson, the topgluon g T , while the additional U(l) symmetry gives rise to a new heavy 
neutral gauge boson, the topcolor Z' T . Topgluons are expected to have a large width, 
so we search for three different widths T = 0.3M, 0.5M, and 0.7M, where M is the 
new particle's mass. For narrow resonances we consider color-octet technirhos from a 
model of walking technicolor ||, a Z' from top color-assisted technicolor [Q], a Z' with 
standard model couplings [Q, and a vector bound state of gluinos (vector gluinonium) 
appearing in supersymmetry models Q. We search for these new phenomena in the 
bb mass spectrum in pp collisions at a center of mass energy \fs = 1.8 TeV. 

A detailed description of the Collider Detector at Fermilab (CDF) can be found 
elsewhere 0. We use a coordinate system with the z axis along the proton beam, 
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transverse coordinate perpendicular to the beam, azimuthal angle 0, polar angle 6, 
and pseudorapidity rj = — lntan(#/2). The silicon vertex detector (SVX), a four- 
layer silicon strip device with radiation-hard electronics, located immediately outside 
the beampipe, provides precise track reconstruction in the transverse plane and is 
used to identify secondary-vertices from b quark decays. The momenta of charged 
particles are measured in the Central Tracking Chamber (CTC), which is inside a 
1.4 T superconducting solenoidal magnet. Outside the CTC, electromagnetic and 
hadronic calorimeters, segmented in rj-cj) towers, cover the pseudorapidity region \rj\ < 
4.2. 

Jets are reconstructed as localized energy depositions in the CDF calorimeters. 
The jet energy E is defined as the scalar sum of the calorimeter tower energies inside 
a cone of radius R = yj (A77) 2 + (A0) 2 = 0.7, centered on the jet direction. The 
jet momentum P is the corresponding vector sum: P = J2 EiUi with being the 
unit vector pointing from the interaction point to the energy deposition Ei in tower 
i inside the same cone. E and P are corrected for calorimeter non-linearities, energy 
lost in uninstrumented regions of the detector and outside the clustering cone, and 
energy gained from the underlying event and multiple pp interactions. Full details of 
jet reconstruction and jet energy corrections at CDF can be found elsewhere 0. 

We define the dijet system as the two jets with the highest transverse momentum 
in an event (leading jets) and define the dijet mass m = J (E~i + E 2 ) 2 — (P~i + i^) 2 - 
The dijet mass resolution is approximately 10% for dijet masses above 150 GeV/c 2 . 
Our data sample was obtained using four triggers that required at least one jet with 
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uncorrected cluster transverse energies of 20, 50, 70 and 100 GeV, respectively. After 
jet energy corrections these trigger samples were used to measure the dijet mass 
spectrum above 150, 217, 292 and 388 GeV/c 2 , respectively At these mass thresholds 
the corresponding trigger efficiencies were greater than 93%. The four data samples 
corresponded to integrated luminosities of 0.087, 2.2, 11 and 87 pb _1 after prescaling. 
We selected events with two or more jets and required that the two leading jets have 
pseudorapidity |^i| < 2 and \r] 2 \ < 2 and a scattering angle in the dijet center-of- 
mass frame | cos 6** | = | tanh[(^i —r] 2 )/2]\ < 2/3. The cos6** requirement ensures 
uniform acceptance as a function of mass and reduces the QCD background which 
peaks at | cos 6 1 * | = 1. To maintain the projective nature of the calorimeter towers, 
the z position of the event vertex was required to be within 60 cm of the center 
of the detector; this cut removed 7% of the events. Backgrounds from cosmic rays, 
beam halo, and detector noise were removed by requiring $t/VJ2 Et < 6 GeV 1//2 and 
E < 2 TeV, where Et is the missing transverse energy, Et is the total transverse 
energy, and E is the total energy in the event. 

To identify jets originating with a 6-quark we require that tracks reconstructed 
in the CTC and SVX form a secondary-vertex, displaced from the event vertex. In 
the secondary-vertex algorithm, described elsewhere 0, we tightened the transverse 
momentum (pt) requirements on the tracks as a function of dijet mass to reduce 
backgrounds at high mass ||. The highest pt track in a reconstructed vertex was 
required to have apr of at least 2 GeV/c for dijets with a mass less than 321 GeV/c 2 . 
This cut was tightened incrementally as a function of dijet mass to 5 GeV/c for dijets 
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with a mass greater than 470 GeV/c 2 . For secondary-vertices formed from only two 
tracks, we additionally required that the tracks have a minimum pt of at least 1 
GeV/c for dijets with mass less than 321 GeV/c 2 , rising to 2 GeV/c for dijets with 
mass greater than 388 GeV/c 2 . We required both jets in the dijet to have a displaced 
secondary-vertex (a b- tag). The efficiency of our 6-tagging algorithm was determined 
from a Monte Carlo simulation of detector response to the decay of a heavy object to 
bb. The simulation was tuned to reproduce the observed tracking efficiency ||. The 
efficiency for 6-tagging a heavy object decaying to bb decreases from 11% to 2.5% 
as the dijet mass increases from 200 to 650 GeV/c 2 . The efficiency decreases as the 
dijet mass increases because the secondary-vertex resolution and acceptance degrade 
as the distance from the event vertex increases and because the tracking efficiency 
degrades as the density of tracks within a jet increases. 

In Fig. we present the inclusive dijet mass distribution for untagged and double- 
6-tagged dijets with \rj\ < 2 and | cos#*| < 2/3. The mass distributions have been 
corrected for the trigger, ^-vertex, and 6-tagging inefficiencies previously discussed. 
We plot the differential cross section versus the mean dijet mass, m, in bins of width 
approximately equal to the mass resolution (RMS~ 10%). The 6-tagged data are 
compared to a smooth parameterization and to a QCD prediction. The QCD pre- 
diction is for direct bb production from the PYTHIA Monte Carlo 0, with the b 
quarks decayed by the CLEO Monte Carlo QQ, and includes a simulation of the 
CDF detector. Direct bb production includes the processes qq — > bb and gg — > bb; 
other bb production processes do not contribute significantly. The QCD simulation 
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used CTEQ2L parton distributions (10| and a renormalization scale /i = Pt- 

To search for resonances, we fit the data to the shape of the bb Monte Carlo 
calculation and a new particle resonance. The bb prediction multiplied by a factor 
of 1.7 fits the data well (x 2 /~DF= 0.61), as shown in Fig. |2| which also shows the 
predicted line shape for narrow resonances and topgluons. Narrow resonances were 
modeled with a PYTHIA simulation of a Z' decaying to bb followed by the CDF 
detector simulation. The mass resolution is dominated by a Gaussian distribution 
from jet energy resolution and a long tail towards low mass from QCD radiation. 
Since the natural width of the Z' is significantly smaller than the reconstructed width, 
these mass resonance curves were used to model the shape of all narrow resonances 
decaying to bb. Topgluons were modeled using a PYTHIA simulation, in which we 
inserted the parton-level sub-process cross section for topgluons 0, followed by the 
CDF detector simulation. We simulated topgluons with width T = 0.3M, 0.5M, and 
0.7M. There is no evidence in the data for either narrow or wide resonances decaying 
to bb. 

Systematic uncertainties on the cross section for observing a new particle in the 
CDF detector are shown in Fig. |[ Each systematic uncertainty on the fitted sig- 
nal cross section was determined by varying the source of uncertainty by ±lcr and 
refitting. The sources of uncertainty presented in Fig. |] are the jet energy scale un- 
certainty, the 6-tagging efficiency, the effect of QCD radiation on the mass resonance 
line shape, the shape of the bb background predicted by QCD, and other sources 
including trigger efficiency, jet energy resolution for narrow resonances, relative jet 



energy corrections between different parts of the CDF calorimeter, and luminosity. 
The dominant systematic uncertainty at low mass is a 3% uncertainty in the jet energy 
scale. The dominant systematic uncertainty at high mass was the uncertainty in the 
6-tagging efficiency, which varied from 14% at low mass to 24% at high mass. The un- 
certainty on the shape of the bb background, from using various parton distributions, 
produced a relatively small uncertainty in the cross section limits compared to the 
dominant sources discussed above, which primarily affect the signal. Other possible 
sources of background, such as mistags or charm, may contribute to the normalization 
but are expected to have a similar shape to bb, and do not contribute significantly 
to the systematic uncertainty. The total systematic uncertainty was found by adding 
the individual sources in quadrature. 

In the absence of evidence for new physics we proceeded to set upper limits on 
the cross section for new particles. For each value of new particle mass in 50 GeV/c 2 
steps from 200 to 750 GeV/c 2 , we performed a binned maximum likelihood fit of the 
data to the background shape and the mass resonance shape. We convoluted each of 
the Poisson likelihood distributions with the corresponding total Gaussian systematic 
uncertainty, and found the 95% confidence level (CL) upper limit presented in Table 
I. 

In Fig. [| we plot our measured upper limit on the cross section times branching 
ratio as a function of new particle mass in 50 GeV/c 2 steps. The limit is compared to 
lowest-order theoretical predictions for the cross section times branching ratio for new 
particles decaying to bb. New-particle decay angular distributions are included in the 
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calculations, and we required \r)\ < 2 and | cos#*| < 2/3 for all predictions. For narrow 
resonances we exclude the color octet technirho in the mass interval 350 < M < 440 
GeV/c 2 . In addition, our limits for narrow resonances are applicable to any particle 
decaying to bb with a width significantly less than our detector resolution of 10%. 
Constructive interference between topgluons and normal gluons causes the signal 
cross section to rise at low bb mass, m, and results in a total cross section integrated 
over all m that is not well defined. To avoid this, the total cross section for a topgluon 
of mass M is defined as the cross section in the region 0.5M < m < 1.5M, for both 
the experimental upper limit and the theoretical prediction. We exclude topgluons 
of width T = 0.3M in the mass range 280 < M < 670 GeV/c 2 , of width T = 0.5M in 
the mass range 340 < M < 640 GeV/c 2 , and of width T = 0.7M in the mass range 
375 < M < 560 GeV/c 2 . 

In conclusion, the measured bb mass spectrum does not contain evidence for a 
mass peak from a new particle resonance. We have presented model independent 
limits on the cross section for a narrow resonance, and set specific mass limits on 
narrow color-octet technirhos and topgluons of various widths. 
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the Republic of China; and the A. P. Sloan Foundation. 
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Table 1: The 95% CL upper limit on the cross section times branching ratio for new 
particles decaying to bb as a function of new particle mass for narrow resonances and 
for topgluons of three different widths (see text). 
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Figure 1: The dijet mass distribution (circles) and double-6-tagged dijet mass 
spectrum (boxes) compared to a fit to a smooth parameterization (dashed curves). 
Also shown is a QCD prediction for bb production (solid curve). 
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Figure 2: The fractional difference between the double-6-tagged dijet mass distribu- 
tion (points) and a background prediction (solid line) is compared to a simulation of 
a 600 GeV/c 2 narrow resonance (dashed curve) and a 600 GeV/c 2 topgluon of width 
T = 0.5M (dotted curve) in the CDF detector. The background has been normalized 
to fit the data and the resonances have each been normalized to the 95% CL upper 
limit on the cross section for a 600 GeV/c 2 resonance. 
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Figure 3: Systematic uncertainties on the cross section upper limit for a) narrow 
resonances and b) topgluons of width T = 0.5M as a function of mass. 



16 




a) Narrow Resonances 
• CDF 95% CL Limit 

Technirho 

Topcolor Z 

Standard Z 1 

Vector 
Gluinonium 



/- 



600 , 
New Particle Mass (GeV/c ) 



200 



b) Topgluons 

• CDF 95% CL Limit 




r/M=o.3 



Excluded: 280 <M< 670 



400 600 

M(g T ) (GeV/c 2 ) 



10 



t 10 

H 
OiD 

u 
PQ 

D 10 



200 



c) Topgluons 

• CDF 95% CL Limit 



r/M=0.5 




Excluded: 340 < M < 640 



400 600 

M(g T ) (GeV/c 2 ) 



200 



d) Topgluons 

• CDF 95% CL Limit' 



r/M=0.7 




400 600 

M(g T ) (GeV/c 2 ) 



Figure 4: The 95% CL upper limit on the cross section times branching ratio (points) 
for a) narrow resonances, and topgluons of width b) T = 0.3M, c) T = 0.5M, and d) 
r = 0.7M is compared to theoretical predictions (curves). 
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